SCIENTIFIC 

REPORTS 



s 




SUBJECT AREAS: 
BIOPHYSICS 
CANCER 
ENGINEERING 
BIOTECHNOLOGY 



Received 
1 August 201 2 

Accepted 
7 January 201 3 

Published 
26 April 201 3 



Correspondence and 
requests for materials 
should be addressed to 
D.W. (wirtz@jhu.edu) 

* A comprehensive list 
of authors and 
affiliations appear at 
the end of the paper. 



A physical sciences network 
characterization of non-tumorigenic and 
metastatic cells 

The Physical Sciences - Oncology Centers Network* 

To investigate the transition from non-cancerous to metastatic from a physical sciences perspective, the 
Physical Sciences-Oncology Centers (PS-OC) Network performed molecular and biophysical comparative 
studies of the non-tumorigenic MCF-lOA and metastatic MDA-MB-231 breast epithelial cell lines, 
commonly used as models of cancer metastasis. Experiments were performed in 20 laboratories from 12 
PS-OCs. Each laboratory was supplied with identical aliquots and common reagents and culture protocols. 
Analyses of these measurements revealed dramatic differences in their mechanics, migration, adhesion, 
oxygen response, and proteomic profiles. Model-based multi-omics approaches identified key differences 
between these cells' regulatory networks involved in morphology and survival. These results provide a 
multifaceted description of cellular parameters of two widely used cell lines and demonstrate the value of the 
PS-OC Network approach for integration of diverse experimental observations to elucidate the pheno types 
associated with cancer metastasis. 

The conversion from a non-tumorigenic state to a metastatic one is of critical interest in cancer cell biology, as 
most deaths from cancer occur due to metastasis \ Typically, we think of the activation of metastasis as one of 
the hallmarks of cancer^ and as a highly regulated, multistep process defined by a loss of cell adhesion due to 
reduced expression of cell adhesion molecules such as E-cadherin, degradation of the extracellular matrix (ECM), 
conversion to a motile phenotype, vascular infiltration, exit and colonization to a new organ site (i.e., intra- and 
extravasation), dormancy, and re-activation. From a physical sciences perspective, metastasis can be viewed as a 
"phase" transition, albeit occuring far from thermodynamic equilibrium^. Though this transition has been the 
focus of much cancer biology research, there is still an incomplete understanding of this phase change, in 
particular, the physical biology of the metastatic state of a cell compared to its pre-malignant state. 
Understanding the physical forces that metastatic cells experience and overcome in their microenvironment 
may improve our ability to target this key step in tumor progression. 

The newly formed Physical Sciences-Oncology Centers (PS-OC) Network, sponsored by and under the 
auspices of the Office of Physical Sciences-Oncology at the National Cancer Institute (OPSO/NCI), is a multi- 
disciplinary network of twelve research centers across the US formed, in part, to test the fundamental hypothesis 
that physical processes (e.g., mechanics, dynamics) play a critical role in cancer initiation and metastasis. The PS- 
OC Network brings analytic techniques and perspectives from the physical sciences to the interpretation of 
biological data and consists of physicists, engineers, mathematicians, chemists, cancer biologists, and computa- 
tional scientists. The goal of the PS-OC Network is to better understand the physical and chemical forces that 
shape and govern the emergence and behavior of cancer at all length scales. The study described in this manu- 
script focused on physical changes associated with metastasis. A controlled set of comparative studies of two cell 
lines that are extensively used as cell models of cancer metastasis and straddle the metastatic transition was 
undertaken by the PS-OC Network. 

The cell lines analyzed were the immortalized human breast epithelial cell line MCF-lOA, representing a non- 
tumorigenic state, and the human metastatic breast cell line MDA-MB-231, representing a malignant state. 
Distinguishing features of the adherent, non-transformed, MCF-lOA cells are their lack of tumorigenicity in 
nude mice, lack of anchorage-independent growth, and dependence on growth factors^. In contrast, MDA-MB- 
231 cells^ form highly malignant, invasive tumors in vivo, are resistant to chemotherapy drugs such as paclitaxel, 
exhibit anchorage-independent growth, and grow independently of growth factors. Although MCF-lOA cells 
have wild-type p53 and MDA-MB-231 cells have mutant p53, both cell lines are negative for the estrogen receptor 
(ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2)^'^. 

To ensure that data generated across the multiple PS-OC laboratories could be integrated, culture guidelines, 
common culture reagents, and the two fully characterized, karyotyped cell lines were distributed to PS-OC 
laboratories. This minimized phenotypic and genotypic drift. After demonstration of growth uniformity, the 
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cells were evaluated by a battery of physical measurements, as out- 
lined in Table 1, encompassing complementary physical, biochem- 
ical, and molecular assays, to establish a metastatic signature across 
multiple length scales, including the molecular, subcellular, cellular, 
and tumor length scales. Novel biophysical techniques interrogated 
classic phenotypic 'hallmark' properties of the two cell lines {e.g., 
morphology, motility, stress responses) and physical cell properties 
(e.g., shear rheology). A novel model-based regulatory network 
approach was used to generate hypotheses of linkages between 
molecular and physical signatures of the cell lines. By interrogating 
this one-of-a-kind dataset, this pilot study provides insight into 
intrinsic differences in the physical properties of metastatic cancer 
cells vs. their non-tumorigenic counterparts, while demonstrating 
the importance of the technologies employed from the physical 
sciences and the value of a network approach to the study of cancer 
biology. 

Results 

In order to generate integrated data across the PS-OC Network, 
the MCF-lOA and MDA-MB-231 cell lines were characterized, 
expanded, and distributed to each PS-OC with common protocols 
and reagents to standardize culturing procedures for each cell line 
(Suppl. Fig. la). Each PS-OC expanded the cells and was required to 
submit predetermined annotations that included phase-contrast 
microscopic images at specified passage numbers, seeding densities, 
and culturing times to ensure phenotypic uniformity of starting 
material (Suppl. Fig. lb). Each of the laboratories then initiated 
experiments using different methodologies to explore cell morpho- 
logy, motility and mechanics, stress responses and survival (drug and 
hypoxia), and molecular networks. 

Comparative cancer cell morphology. Cancerous cells are often 
described as having an altered appearance and morphology; cancer 
pathologists and oncologists routinely use cell and nuclear morph- 
ology to stage cancer and propose treatments. Clinical studies have 
linked properties of tumor cell mass and patient survival to variations 
in individual cells^ and nuclear morphology^. In particular, breast 
cancer cell line morphologies have been correlated to invasiveness 
and gene expression profile In this study, morphological 



parameters of the MCF-lOA and MDA-MB-231 cells, including 
cell and nuclear shape, nuclear disorder strength, membrane lipid 
raft abundance, and cell growth in response to matrix stiffness were 
characterized. To determine the two- (2D) and three-dimensional 
(3D) properties of the two cell lines, MCF-lOA and MDA-MB-231 
cells were grown as monolayers and imaged in 2D by differential 
interference contrast microscopy (Fig. la). The MDA-MB-231 cells 
had a higher width-to-length aspect ratio of 1 : 1 .77 compared with the 
smaller, rounder MCF-lOA cells' width-to-length aspect ratio of 
1:1.39 (Fig. la). Three-dimensional cellular and nuclear shapes 
were determined using single- cell optical computed tomography 
and volumetric image analysis^\ Representative pseudo-color 
renderings of 3D cell volumes are shown in Fig. lb. Detached 
MDA-MB-231 cells exhibited consistent kidney-bean shaped 
nuclei with a relatively consistent extent of concavity, whereas 
nuclei of MCF-lOA cells were more heterogeneous in shape. 
Nuclear sphericity parameters calculated from phase- contrast 
images were 1.42 and 1.39 for the MDA-MB-231 and MCF-lOA 
cells, respectively, demonstrating the significantly non-spherical 
shapes of the nuclei of both cell types (Fig. lb iii). This is contrary 
to the popular notion of suspended cells having spherical nuclei and 
is a novel observation of nuclear shape for these two otherwise well- 
characterized cell lines. 

Partial wave spectroscopic (PWS) microscopy is a recently intro- 
duced high- resolution approach to characterize cancer cell morpho- 
logy and the nanoscale architecture of the nucleus PWS employs 
back- scattered light to quantify the spatial variations of the refractive 
index or macromolecular mass density. Measurements of macromo- 
lecular compaction are reported as the nuclear disorder strength, L^, 
where an increase in the disorder strength of the cell nanoarchitec- 
ture is used as an indicator of early events in carcinogenesis^^. 
Nuclear regions of the MDA-MB-231 and MCF-lOA cell lines were 
analyzed and calculated for each cell line and normalized to MCF- 
lOA values as baseline, as illustrated in Fig. Ic (PWS images are 
shown on the right and the bright field reflectance (BFR) images, 
from which PWS images were computed, are on the left). Metastatic 
MDA-MB-231 cells exhibited a statistically significant (p < 0.0004) 
80% increase in nuclear disorder strength relative to the non-tumori- 
genic MCF-lOA cells. This difference in nuclear indicates that the 



Table 1 Cell physical parameters, methods and measurements 


Cell physical parameter 


Technology name/ Physical method 


Measurement 


Morphology 


Cell and nuclear shape, volume 

Nuclear architecture 

Cell growth as a function of matrix 

stiffness 
Cell surface 


2D: Differential interference contrast microscopy 
3D: Optical computed tomography 

Partial wave spectroscopy 
Immunofluorescence; confocal microscopy 

Total internal reflection fluorescence; epifluorescence 
microscopy 


2D: Shape; length 

3D: Volume; characteristic nuclear shape (nuclear 

sphericity parameter) 
Nuclear disorder strength (L^) 
Cell proliferation; cell morphology 

CD44 expression patterns; lipid raft distribution 


Motility and Mechanics 


Cell motility 

Cell adhesion and rolling 
Mechanical flexibility 
Internal fluidity 

Endogenous force generation 


1 D, 2D, and 3D motility assays 

HA Micropatterns; flow chamber adhesion assay 

Atomic force microscope-based nano-indentation 

Ballistic injection nanorheology 

Traction force microscopy 


Speed; radial displacement 

Cell binding to HA pattern; rolling velocity 

Elastic modulus 

Mean square displacement 

Tension maps, force magnitudes 


Stress Response and Survival 


Hypoxic conditions 
Chemical stress 


Microscopic imaging of 3D cultures (intracellular Ca^^ 

and DNA) 
As above 


2D Cell viability; 3D O2 consumption; CEA 

expression 
Cell viability 


Abbreviations: 2D: 2-dimensional; 3D: 3-dimensional; CEA, carcinoembryonic antigen; HA, hyaluronic acid. 
Footnote: Additional PS-OC methods are noted in the results and given in the supplementary information. 
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Figure 1 | Comparative cell morphology, (a) Differential interference contrast (DIG) microscopy, (i, ii) Left: Volume rendering from DIG micrographs 
of each cell type (gray, H&E stained). Right: Same as left with EFM images of DAPI stained nuclei (blue) superposed, (iii) Aspect ratios of cell bodies 
(mean ± s.e.m.). (b) 3D cytometry. (i,ii) Pseudo-colored volume rendering of suspended and fixed H&E stained cells imaged by optical cell GT. 
Gytoplasm is grey and nucleus is blue, (iii) Nuclear sphericity of the cell nuclei (mean ± s.e.m.). (c) Nuclear disorder strength, (i, ii) Left: Bright field 
reflectance (BFR) images. Right: PWS microscopic images. Golor shows the magnitude of the nuclear disorder strength (L^) (low: blue, high: red). L^ 
values normalized to 1.0 for MGF-lOA cells, (iii) Nuclear disorder strength (mean ± s.e.m.). (d) Substrate stiffness, (i, ii) Gonfocal immunofluorescence 
of cells grown for 15 days in 3D on soft (Left, 75 Pa) and hard (Right, 6000 Pa) reconstituted basement membrane-conjugated polyacrylamide gel 
matrix. Gells stained for Ki-67 cell proliferation marker (red) and DNA using DAPI. (iii) Fraction of Ki-67 positive cells as function of substrate stiffness 
(mean ± s.e.m.). (e) GD44 distribution, (i, ii) GD44 distribution visualized by anti-GD44 antibodies using (Left) epifluorescence microscopy (EFM) and 
(Right) total internal reflection fluorescence (TIRF). (iii) Fluorescent area in \im^ calculated from TIRF images (mean ± s.e.m.). Bar graph compares 
significant difference in immunofluorescence intensity between TIRF images for MDA and MGF cells, not the EFM images, (f) Lipid raft distribution, 
(i, ii) Left: Lipid raft distribution visualized by anti-GT-B antibodies using EFM. Right: raft distribution visualized by anti-GT-B antibodies using TIRF. 
(iii) Fluorescent area in [im^ calculated from TIRF images (mean ± s.e.m.). All scale bars are 5 |am. All p- values are indicated according to the Michelin 
guide scale (p < 0.001: [***]; 0.001 <p< 0.01: [**]; 0.01 <p< 0.05: [*]; 0.05 < p: ns). 



metastatic and non-tumorigenic cells have very different nuclear 
architectures. 

There is increasing evidence of a functional relationship between 
tissue rigidity and tumor progression; indeed, tumors are often stiffer 
than normal tissues^^. An investigation of whether cell matrix stiff- 
ness modulates cell growth and morphology was carried out with the 
two cell lines using ECM-crosslinked polyacrylamide gels of varying 
stiffness (Fig. Id). The substrate stiffness was measured by shear 
rheology (the study of the flow of matter) When MCF- IDA cells 
were grown for more than 15 days on a soft matrix (with sub- 
strate stiffness of either 75 Pa or 140 Pa), the cells formed well- 
differentiated, hollow, growth- arrested acini. When the cell matrix 
stiffness was increased to 6000 Pa, the stiffer substrate induced an 
increase in cell growth but cells formed poorly differentiated mono- 
layers (Fig. Id) with nearly 70% of cells highly proliferative as 
evidenced by Ki-67 positive- staining. Conversely, MDA-MB-231 
cells exhibited similar morphology and growth rate regardless of 
substrate, suggesting that these cells have a significantly reduced 
mechano- sensitivity. This insensitivity to matrix stiffness may be 



beneficial during metastasis as metastatic cells encounter different 
substrates during the dissemination process, and the lack of sensitiv- 
ity to matrix stiffness allows these cells to proliferate in a wide variety 
of environments. 

The expression and distribution of cell surface glycoproteins are a 
major determinant of cellular response to the microenvironment^^. 
One of these glycoproteins, CD44, is a hyaluronic acid (HA) receptor 
involved in cell-cell adhesion and cell-matrix interactions, which are 
critical to the metastatic process. Variant isoforms have been assoc- 
iated with cancer metastasis, particularly in tumors originating from 
epithelia, including breast cancer^^'^^. Investigation of the surface 
receptor expression on MCF-lOA and MDA-MB-231 cells using 
epifluorescence microscopy (EFM) images (Fig. le) revealed that 
the macroscopic CD44 distributions of the two cell lines did not 
differ significantly (p > 0.75) (data not shown). This confirms a 
previous report using confocal microscopy, which showed a similar 
localization of CD44 on MDA-MB-231 cells^^ In contrast, total 
internal reflection fluorescence (TIRF) images of the same cells 
(Fig. le) indicated that CD44 was abundantly present at the point 
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of contact between the cell membrane and the substrate surface in the 
MCF-lOA but was observed at markedly lower levels (by a factor of 
~4.5) in the MDA-MB-231 cells (Fig. le iii). Interestingly, CD44 
expression and distribution were not affected in either MCF-lOA 
or MDA-MB-231 cells when cells were grown in 3D on substrates 
of differing stiffness. CD44 was mainly present at cell- cell junctions 
and was evenly distributed along the plasma membrane in both cell 
types (data not shown). When the cells were treated with a fluor- 
escent probe specific to membrane lipid rafts, enhanced labeling of 
the MDA-MB-231 cells relative to MCF-lOA cells was observed 
(Fig. If), suggesting the presence of lipid rafts only in the metastatic 
breast cancer cells. This finding is consistent with reports of elevated 
levels of cholesterol-rich lipid rafts in breast and prostate cancer cells 
compared to their normal counterparts^^. Taken together, these 
experiments suggest that although the general topography of the 
two cell lines are similar, the surface presentation of the cell adhesion 
protein CD44 and formation of lipid raft domains containing addi- 
tional cell surface receptors are significantly different. 

Comparative cancer cell motility and mechanics. Cell migration. 
One of the distinguishing hallmarks of metastatic cells is their 
capacity to steer through multiple physical microenvironments 
such as the ECM of the stromal space and, following intravasation, 
along vasculature walls. Traditionally cell motility studies have been 
performed in 2D environments (i.e., flat substrates). Here, cell 
motility was evaluated in one-dimensional (ID), 2D, and 3D 
environments^^"^\ These multidimensional experiments revealed 
that regulation of cell speed and maximum displacement were 
critically dependent on the dimensionality of the environment. In 
a ID environment where a cell is constrained to move only forward 
or backward (13 |im wide by 25 |im deep silicon -etched fibronectin- 
coated channels), MCF-lOA cells traveled more than three times 
faster than MDA-MB-231 cells (Fig. 2a iii). However, consistent 
with their metastatic potential, MDA-MB-231 cells traveled farther 
(based on radial displacement from their original position) than 
MCF-lOA cells along the same ID channels (Fig. 2a iii). In a 2D 
environment, MCF-lOA cells moved in a circular or pin- wheel 
style of motility — the leading edge swung in an arc while the 
lagging edge often remained pinned in place (data not shown). 
MDA-MB-231 cells moved more linearly, though more slowly, and 
MCF-lOA cells were found to travel farther than MDA-MB-231 cells 
on the 2D collagen matrix (Fig. 2a iii). MCF-lOA cells also exhibited 
slightly faster migration in a wound healing assay (see Suppl. Fig. 5). 
MCF-lOA cells also moved in a circular motion around the 3D void 
in which they were embedded, but MD A-MB-23 1 cells again traveled 
farther than MCF-lOA cells when embedded inside 3D collagen 
matrices (data not shown). Thus, given the dimensional con- 
straints of the cellular environment, the non-tumorigenic cells 
tended to move faster than the metastatic cells, but remained 
within a limited circular area, whereas MDA-MB-231 cell motility 
was linear and did not exhibit the same distance limitations. 

ECM matrix components are known to potentiate breast tumor 
metastasis by enhancing cell invasion^^, therefore, the functional 
relationship between cell motility and ECM molecules such as lami- 
nin (a primary ECM component of the basement membrane in 
breast tissue), hyaluronic acid (HA, an anionic non-sulfated glyco- 
saminoglycan spatially distributed in the ECM), and cell surface 
adhesion molecule E-selectin were investigated. Traction force 
microscopy (TFM) was used to measure the pulling forces exerted 
by single cells on substrates with different concentrations of laminin. 
MDA-MB-231 cells generated more traction forces than MCF-lOA 
cells at all tested laminin concentrations (0.1 |ig/ml, 10 |ig/ml, and 
50 |ig/ml) (Fig. 2b iii, intermediate concentration data not shown), 
suggesting that laminin interactions may play a role in promoting 
breast tumor cell aggressiveness. Cell adhesion to HA substrates and 
CD44 expression in the two cell types were compared by growing 



cells on micropatterned arrays with covalently linked HA (Fig. 2c). 
Immunofluorescence and flow cytometry experiments confirmed 
that MDA-MB-231 cells had approximately two-fold more 
CD44 expression than MCF-lOA cells (Fig. 2e iii). Nevertheless, 
MCF-lOA cells adhered preferentially to HA, whereas MDA-MB- 
231 cells showed no preference (Fig. 2c iii). The lack of HA adhesion 
preference of the metastatic cells is likely due to the lack of CD44 in 
the membrane. This is evident in the TIRF data (Fig. le iii). Whereas 
the EFM analysis (data not shown) shows 1:1 bulk CD44 levels, by 
TIRF, levels of CD44 are 2:1, a difference presumably due to the 
presence/absence of external HA in the environment. 

The comparative role of E-selectin in the cell migration of meta- 
static and non-malignant cells may play a critical role in the ability to 
adhere to vasculature walls for efficient extravasation to secondary 
organs^^. For example, in trans -endothelial migration of circulating 
tumor cells (CTC), endothelial cell surface ligands increase the 
adhesive forces and hence the residence time and extravasation of 
CTCs from the circulatory system. Fluid shear forces due to blood 
flow can counteract these forces. Shear stresses from 1-4 dyn/cm^ are 
typical of veins, whereas higher stresses ranging from 4-30 dyn/cm^ 
can occur in arteries\ Here, flow-chamber experiments showed that 
MDA-MB-231 cells neither adhered nor rolled on either 5 |ig/mm^ 
or 10 |ig/mm^ E-selectin-coated surfaces over a range of physio- 
logical wall shear stresses from 1-8 dyn/cm^ (data not shown). In 
contrast, MCF-lOA cells adhered and rolled on 5 |ig/mm^ E-selectin- 
coated surfaces over the entire range of stresses ranging from 1 to 
8 dyn/cm^ (Fig. 2d iii). These results differ from those of Zen et aU^, 
who found that MDA-MB-231 cells were able to traverse a model 
endothelial monolayer in the absence of shear forces in a CD44- and 
E-selectin- dependent manner^^. Given the expression of CD44 on 
MDA-MB-231 cells, the CD44/E-selectin binding interaction (affin- 
ity and/or expression) on the cell surface may be insufficient to 
increase the residence time in the presence of significant shear forces, 
and other ECM ligands or adhesion molecules may compensate in 
vivo. In spite of an apparent lack of HA and E-selectin binding by 
MDA-MB-231 cells, both MDA-MB-231 and MCF-lOA cells exhib- 
ited an ECM deposition similar to that of Nuff fibroblast cells (data 
not shown), however the ECM structures appear to be distinct 
(Fig. 2e i). 

Cell mechanics. Consistent with observations of amoeboid move- 
ment of invasive cancer cells, the ability of a cell to move through 
multiple tissue compartments, often via small portals, relies on 
amoeba-like deformability. The mechanical deformability of MCF- 
lOA and MDA-MB-231 cells was measured by an atomic force 
microscope (AFM) aligned with a confocal fluorescence microscope 
lens (AFM-CLSM) for fluorescence lifetime imaging measurements 
(FLIM) (Fig. 2f)^^. Cells stained with nuclear and nucleolar dyes were 
indented at distinct points over the lamella, nucleus, and nucleoli. 
The resulting force-indentation curves were fitted to a modified 
Hertz model in 100-nm intervals, yielding depth- dependent elastic 
moduli. At shallow indentation depths, the two cell lines had similar 
elastic moduli (—200 Pa). Cytoplasmic and nuclear stiffness of 
MDA-MB-231 cells increased only slightly with increasing indenta- 
tion depth, whereas cytoplasmic and nuclear stiffness of MCF-lOA 
cells both increased - though with different strain-hardening profiles 
- about four-fold to an elastic modulus of — 1.6 kPa. At indentation 
points over nucleoli, both cell lines showed stiffening with increasing 
depth, with MDA-MB-231 and MCF-lOA cells rising to elastic mod- 
uli of — 1 kPa and —1.5 kPa, respectively, at 0.8 |im. These findings 
suggest that mechanical loads are transduced through the cytoskele- 
ton differently in the two cell lines. The decreased elastic modulus 
and increased deformability of MDA-MB-231 cells is consistent with 
their ability to traverse narrow matrices. 

Cytoplasmic viscoelasticity is also an important feature of amoeb- 
oid movement. Subcellular viscoelasticity was investigated by bal- 
listic injection nanorheology (BIN; movement of nanoparticles in a 
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Figure 2 | Motility and mechanics, (a) Three substrates used to monitor cell motility: 2D collagen-coated glass substrate, 3D collagen matrix, and ID 
fibronectin-coated microchannels (13 \im wide, 25 [im deep) etched in silicon. Motion tracking based on time-lapse imaging, (i) Schematic, (ii) Tracking 
data, (iii) Speed (|am/min) and maximum invasion distance (mean ± s.e.m.). (b) TFM quantified traction stresses exerted by cells on 5kPa 2D 
polyacrylamide substrate mimicking mammary tumor stiffness by measuring displacement of embedded fluorescent polystyrene beads. Phase image 
followed by fluorescent images of bead field under stressed and unstressed (post-trypsinization) conditions. Bead displacement yields magnitude ( ITI ) 
and distribution of traction stresses^^. (i) Schematic, (ii) Tension maps (Left); phase images (Right); MCF-lOA (Top); MDA-MB-231 (Bottom). Scale 
bars: 50 |am. (iii) Force magnitudes exerted by cells at different surface laminin concentrations (mean ± s.e.m.). (c) Microprinted Covalent HA array, (i) 
Schematic, (ii) CD44 expressing MCF-lOA and MDA-MB-231 cells (red) attached to FL-HA micro-patterned substrates (green) after 24h culture. Scale 
bars: 50 |am. (iii) Cells attached to HA squares (mean ± s.e.m.). (d) Cells rolling on E-selectin surfaces, (i) Schematic, (ii) Phase images. Scale bars: 50 \im. 
(iii) Rolling velocities and numbers of MCF-lOA cells captured on surface under physiological wall shear stresses (mean ± s.e.m.). (e) (i) SEM reveals 
distinct ECM structures deposited by MCF-lOA (Left; scale bar: 4 |am) compared to MDA-MB-231 (Right; Scale bar: 5 |am). (ii) IF imaging of cells 
stained with fluorescein-tagged HA demonstrate expression of HA in both MCF-lOA and MDA-MB-231. Scale bars: 50 [im. (iii) Flow cytometry shows 
higher CD44 expression in MDA-MB-231 compared to MCF-lOA (left: histogram; right: quantified MFI). (f) AFM probe aligned with confocal 
fluorescence lifetime microscope scans points of interest over cytoplasm, nuclei, and nucleoli. Force-indentation curves used to calculate elastic moduli, 
(i) Schematic, (ii) Curves (middle) and corresponding images (top, bottom). Scale bars: 10 [im (top); 2 [im (middle, horizontal), 0.2nN (middle, 
vertical); 4 |am (bottom), (iii) Depth -dependent elastic moduli (mean ± s.e.m.). (g) Fluorescent nanoparticles injected into cells and trajectory 
monitored over time, (i) Schematic, (ii) Cell monitored in real time. Inset: nanoparticle trajectory. Scale bars: 10 [im (main); 0.2 [im (inset), (iii) MSD 
values over cumulative time (mean ± s.e.m.). All p-values indicated by Michelin guide scale (p < 0.001:[***]; 0.001 <p < 0.01:[**]; 0.01 <p< 0.05:[*]; 
0.05 <p:ns). 
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viscoelastic material)^^'^^. Fluorescently labeled, 100 nm diameter 
microspheres were ballistically injected into the cytoplasm and, after 
overnight incubation, their random intracellular displacement was 
followed over time (Fig. 2g). Mean squared displacement (MSD) of 
the microspheres was obtained from the 20 s trajectory of each 
microsphere with 30 ms temporal resolution. Consistent with the 
above AFM results, the ensemble-averaged MSD of microspheres 
in MDA-MB-231 cells was greater than that of MCF-lOA cells, indi- 
cating that the MDA-MB-231 cytoskeleton was substantially softer 
than that of MCF-lOA cells. 

Comparative cancer cell stress response and survival. Preferential 
survival under stressful conditions is a characteristic of metastatic 
cancer cells; therefore, the two model cell lines were subjected to 
external stresses characteristic of the tumor microenvironment, 
such as hypoxia (in 2D and 3D) and low pH. Cell viability, 
recovery, oxygen consumption (of single cells and populations), 
and expression of the surface biomarker carcinoembryonic antigen 
(CEA) were measured. 

Hypoxia. Within the primary tumor, oxygen availability can vary 
dramatically due to location within the tumor and its vascularity; 
oxygen availability also varies temporally {e.g., due to clots). The 
effects of hypoxia on cell viability, growth and recovery, and oxygen 
consumption were tested in a 2D environment as well as in a patho- 
logically relevant 3D culture environment. In a hypoxic (1% O2) 2D 
environment, populations of both cell lines experienced relatively 
small decreases in viability over a span of three days (Fig. 3a, top). 
In a 3D culture environment, viability of both cell lines was appar- 
ently more sensitive to hypoxia (Fig. 3a, bottom). The normalized 
MCF-lOA cell viability under ambient and hypoxic conditions at day 
6 were 100 (± 20) % and 40 (± 10) %, respectively. The normalized 
MDA-MB-231 cell viability under ambient and hypoxic conditions 
at day 6 were 150 (± 30) % and 76 (± 20) %, respectively. The 
apparent differences of viability in hypoxic 2D and 3D environments 
require further study but may involve oxygen permeability (surface 
vs. bulk effects in 3D), dead cells remaining caught in 3D scaffolding, 
or changes in cell adhesion. O2 media concentrations were also mea- 
sured in the 3D cultures and the normalized bulk oxygen consump- 
tion rate at day 6 was determined for both cell lines (Fig. 3b, left). 
Under ambient conditions, MCF-lOA and MDA-MB-231 cells 
exhibited similar O2 consumption rates; however, their response to 
hypoxic (1% O2) treatment differed with no changes in MCF-lOA O2 
consumption, whereas MDA-MB-231 cells dramatically reduced O2 
consumption (four-fold) when cultured in hypoxia (Fig. 3b, left). The 
observed reduction of O2 consumption by MDA-MB-231 cells may 
be due to an increased plasticity or adaptability related to their 
tumorigenic potential. An alternative explanation may be the preex- 
istence of variant MDA-MB-231 cells in the population with reduced 
O2 metabolism that are positively selected for under hypoxic stress. 
To examine these possibilities, oxygen consumption rates of single 
cells isolated in hermetically sealed chambers were measured in nor- 
moxic (17% O2) conditions (Fig. 3b, right). In these measurements, 
MDA-MB-231 cells had a lower mean OCR but comparable hetero- 
geneity (mean: 2.1 fmol/min, CV: 0.72) than the MCF-lOA cells 
(mean: 4.1 fmol/min, CV: 0.82). The histogram suggests a glycolytic 
subpopulation of the metastatic cells may become dominant in hyp- 
oxic conditions, but further studies are required. This hypothesis is 
sketched out in Fig. 3c. In addition, measurements of the expression 
of selected cell surface proteins thought to be involved in metastasis, 
such as CEA (Fig. 3d), PCLP, and CD44 (Suppl. Fig. 2d), as a function 
of a normal or hypoxic environment showed that MDA-MB-231 
cells increased expression of CEA in response to hypoxic conditions. 
Future investigations may clarify the effect of hypoxic conditions on 
the expression of these proteins. 
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Figure 3 | Comparative cell stress responses, (a) Viability under hypoxia in 
2D and 3D. Top: (2D) Cells grown in weUs (triplicate) for 3 days in 1% O2. 
Viability was determined every 24 h and imaged with an inverted microscope. 
Cell viability (mean ± s.e.m.) normalized to day 1 samples. Bottom: (3D) 
DNA content per scaffold (normalized to day 1 samples) over 6 days growth 
in normoxic (17% O2) and hypoxic (1% O2) conditions in 3D culture 
(alginate discs), (b) Oxygen consumption rates. Left: bulk OCR (normalized 
to DNA content, mean ± s.e.m.) of cells after 6 days growth in normoxic 
(17% O2) or hypoxic (1% O2) conditions in 3D culture (alginate discs); 
Right: histogram of OCR of single cells measured in hermetically sealed 
chambers (17% O2). (c) Schematic of a hypothetical model in which 
phenotypic diversity of MDA-MB-231 cells is relatively enhanced with respect 
to MCF-lOA due to enhanced population recovery after hypoxia-induced cell 
death, (d) Carcinoembryonic antigen (CEA). Mean fluorescence intensity 
minus isotype (MFI) of cells grown in 17% or 1% O2 (mean ± s.e.m.). 
(e) pH-induced stress. Percentages of viable cells grown in media with pH 6.8 
(mean ± s.e.m.). (f) Paclitaxel-induced stress. Percentages of viable cells after 
24, 48, and 72 h incubation with various concentrations of paclitaxel (mean 
± s.e.m.). AH p- values are indicated according to the Michelin guide scale 
(p < 0.001: [***]; 0.001 < p < 0.01: [**]; 0.01 < p < 0.05: [*]; 0.05 < p\ ns). 
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Chemical stress. Changes in the interstitial chemistry of the micro- 
environment can also have a profound impact on cancer cell res- 
ponse and metastasis. Low interstitial pH is typical in the tumor 
microenvironment due to increased lactic acid secondary to anaero- 
biosis within anoxic tumors. When cells were grown under acidic 
conditions, pH 6.8 vs. neutral pH 7.4, and cell proliferation was 
monitored at fixed time points (microscopic imaging of calcein- 
AM and ethidium homodimer-1 staining every 24 h for 72 h), 
MCF-lOA viability dropped to 70% of control by 24 h and remained 
there through 72 h (Fig. 3e), whereas viability of the metastatic cell 
line was not significantly affected. Additionally, viability of the cell 
lines was measured after 24 h, 48 h, and 72 h exposure to a range of 
doses of the mitotic inhibitor paclitaxel (Fig. 3f). Although a clear 
dose-dependent decrease in viability was observed in both cell lines, 
the MDA-MB-231 cells were considerably less sensitive to the drug 
than the MCF-lOA cells. 

Molecular network signatures for morphology, motility, and 
stress. The goals of our studies were to broaden the set of features 
that differentiate between metastatic and non-tumorigenic cells to 
include physical/mechanical properties and to generate hypotheses 
about how molecular-scale factors {e.g., protein and transcript level 
changes) impact or are impacted by these properties. The following 
molecular network analysis attempts to identify putative molecular 
origins of the data from Fig. 1 and Fig. 2. Notably, in addition to the 
transcriptomic data, proteomics datasets from each cell line were 
collected both unperturbed and perturbed by various antagonists 
that could potentially alter parameters such as stress response {e.g. 
paclitaxel) and motility and morphology (ROCK Inhibitor Y-27632). 
This allowed us to further refine our estimates of definition of genes 
from the network that were putatively connected to particular 
phenotypic characteristics. To connect the molecular and the 
biophysical, we proceeded in three steps: 1) we identified relevant 
transcriptome data (224 experimental conditions from the Gene 
Expression Omnibus) for the two cell lines and generated 
quantitative differential proteomics data (see Methods); 2) we 
derived a computational model of cellular regulation; and 3) we 
identified subsets of that model that are likely related to physical 
properties of the cells. 

Our transcriptional regulatory network model (Suppl. Fig. 3) con- 
tains 1866 genes, 220 of which encode transcription factors (TFs) and 
highlights factors that differentiate between metastatic and pre- 
malignant cell types (see Methods). In this model, each gene is repre- 
sented as a 'node'. If a gene's abundance is regulated by another gene, 
this is denoted with an 'edge' between those genes. In Fig. 4a, we show 
the transcription factor subset. Yellow nodes or edges indicate spe- 
cificity to MCF-lOA cells, whereas blue nodes or edges are more 
specific to MDA-MB-231 cells. This network shows a substantial 
bias towards MDA-MB-231 specific nodes and edges. 

We interrogated our derived total network (Suppl. Fig. 3) to 
identify sub-networks whose regulation was both highly differen- 
tiated between cell types and likely connected to physical cell prop- 
erties (morphology, motility, and stress). Proteomics measurements 
were overlaid onto these networks. 

In Figs. 4b-d, we show some of the 1-hop sub-networks {i.e., one 
degree of separation between genes) selected from the millions of 1- 
hop sub -networks identified. These particular networks contain one 
or more interesting features: connections are specific to one cell line; 
the genes are of interest to the PS-OC community and/or implicated 
in cell biophysical properties; and/or the genes correspond to differ- 
entiating proteomics data. 

As noted above, there were numerous morphologic differences 
between our cells, including factors like width-to-length ratio and 
nuclear disorder factors. Notably, for the morphology sub-network 
(Fig. 4b), we observe strong connections among FBNl (encoding 
fibrillin, a major extracellular microfibril in connective tissues). 



ZEBl (encoding a zinc finger transcription factor that represses T- 
lymphocyte-specific IL2 gene), and TWIST 1 (encoding a basic helix- 
loop-helix transcription factor). The inferred network connects these 
genes, whose products are known to impact the morphology of cells. 
FBNl is of particular interest because it has a large number of 
MDA-MB-231 -specific edges. It is a connective protein that provides 
structural support for surrounding tissue. The large number of spe- 
cific edges suggests that FBNl plays a role in the deformability of the 
MDA-MB-231 cells. The network analysis further infers that FBNl 
has regulatory interactions with gene products involved in cellular 
differentiation, adhesion, structure, and integrity, such as TGFPS 
(encoding transforming growth factor P3), MMP2 (encoding matrix 
metalloproteinase-2), LOX (encoding lysyl oxidase), and ACTA2 
(encoding smooth muscle aortic a-actin). TGFPS, MMP-2, and 
LOX have previously been shown to be involved in tumor growth 
and metastasis^^'^^, whereas to our knowledge ACTA2 has not. 

An analysis of motility network signatures in MCF-lOA and 
MDA-MB-231 cells identified lTGt^4 (encoding integrin P4) as a 
central node (Fig. 4c). Integrin P4 is a member of the integrin family 
necessary for cell-matrix and cell-cell adhesion and for motility^^'^°. 
Our model suggests that it is a potential regulator of many gene 
products involved in adhesion, migration, and invasion. This par- 
ticular network also contains a large number of differentially abund- 
ant gene products based on proteomics, most notably SERPINB5 
(encoding serpin peptidase inhibitor, clade B, member 5), which is 
upregulated in the non-maUgnant MCF-lOA cells and is known to 
act as a tumor suppressor and block the metastatic properties of 
mammary tumors, further suggesting a role in differential regulation. 

In the survival network, we focused initially on HIFIol (hypoxia- 
inducible factor- la), as it encodes a transcription factor known to be 
involved in cell response to hypoxia^ \ We focused on this gene 
because it is over- expressed in the metastatic cells able to survive 
in areas of low O2 within a tumor (Fig. 4d). Products of genes in this 
regulatory network ultimately affect the survival of a cell through 
proliferation or apoptosis. Our model suggests that c-Met has reg- 
ulatory relationships with both HIFla and ITGP4 in the transcrip- 
tion factor regulatory network, consistent with the reports from 
Schelter et aV^ and Giancotti^^ respectively. In addition, we suggest 
that HlFloL regulates LOX, which has been shown to be critical for 
hypoxia-induced metastasis^^"^^. 

Discussion 

This trans PS-OC laboratory network study demonstrates that cell 
lines can be cultured with genotypic and phenotypic uniformity 
across geographically and technologically disparate institutions to 
be used in comparative studies. The results from each laboratory 
can be pooled with confidence, given the known uniformity in cell 
maintenance and thus the phenotypic and genotypic make-up of the 
cells. This is in contrast to typical comparisons with results from 
literature that are fraught with the pitfalls inherent in the use of 
different experimental conditions and materials. In the present work, 
the physical differences between a non-tumorigenic cell line (MCF- 
lOA) and a metastatic cell line (MDA-MB-231) that are both exten- 
sively used in research were characterized. The integrated results 
generate a composite and multifaceted picture (detailed in Table 2) 
that is more complete than one or a few experimental approaches 
would have provided. In proposing to develop a laboratory network 
approach to study the physical properties of cancer cells, there was 
concern about the inherent mutator genotype of the cell lines as well 
as phenotypic plasticity due to variant culturing conditions, which 
could result in differing cell lineages and phenotypes potentially 
confounding interpretation and integration of results across the 
twelve PS-OCs. Consequently, considerable work was done to ensure 
that the two cell lines, culture media, and growth conditions were 
standardized and that cell growth and overall morphology were 
documented at regular intervals for quality control. Subsequently, 
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Figure 4 | Comparative molecular signatures for morphology, motility, and stress, (a) The largest connected subnetwork of transcription factors from 
the master network (Suppl. Fig. 3) with nodes colored to provide a "summary" of the entire network. Node size shows the number of edges (connecting 
lines) in the master network that were above a cutoff for specificity to either cell line. Larger nodes have more cell-line-specific edges; the largest, IKZFl, 
has 67 edges above the threshold. Node color is determined by the ratio of above-cutoff edges specific to MCF-lOA vs. MDA-MB-231, with yellow 
denoting more MCF-lOA edges and blue more MDA-MB-231 edges. Nodes with many edges specific to one cell line or the other are therefore large and 
brightly colored, such as IKZFl or COPS2. (b-d) One-hop networks from transcription factor regulators (A) to their targets ( O). Each gene is 
represented as a 'node'. If a gene's abundance is regulated by another gene, this is denoted with an 'edge' between those genes. Color of an edge indicates 
the specificity of that regulatory relationship to either MCF-lOA cells (yellow) or MDA-MB-23 1 cells (blue). Relationships that are equally present in both 
cell types are demarked grey. Node border color indicates differential proteomics results. Yellow border nodes are upregulated in MCF-lOA cells. Blue 
border nodes are upregulated in MDA-MB-231 cells. Grey bordered nodes were quantified and found to be equivalent in both cell types, (b) Morphology 
network. The 1-hop morphology network from FBNl and TWISTU LOX and LOXLU both putatively regulated by FBNl. Both FBNl and TWISTl are 
putatively regulated by ZEBl. Also shown are the large number of MDA-MB-231 edges from FBNl and a fairly even distribution of edges from ZEBl. 
(c) Motility network. The 1-hop network from ITGB4. ITGB4 is itself a gene of interest and is inferred to regulate EGFR and several laminins. (d) Stress 
response network. The 1 -hop network from HIFIA, a transcription factor and gene of interest. It is putatively regulated by MET (upper triangle) , which is 
also inferred to regulate ITGB4. HIFIA putatively regulates two more genes of interest, LOX (also a putative target of FBNl and SATB2). 

each lab employed their own unique procedures for cell character- phenotypic plasticity of the metastatic MDA-MB-231 cells compared 

ization as provided in Table 1 and the supplemental information. with the non-tumorigenic MCF-lOA cells is illustrated at the 20 nm 

Integration of the results across the study highlights the similar- scale, at which an increase in nuclear disorder strength indicative 

ities and differences between the cell types. For example, the apparent of more variable nuclear heterochromatin structure was seen in 
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MDA-MB-231 cells, and at the jim scale, at which the cellular 
deformability of the metastatic cells was found to be four-fold more 
elastic than the non-tumorigenic cells as measured by AFM 
(Table 2). These results are consistent with the role of the cytoskele- 
ton in malignancy. The cytoskeleton functions not only as the cell 
backbone but also plays an active role in cell division, differentiation 
and signal transduction; modifications in the cytoskeleton have been 
associated with early carcinogenesis^^. In addition, the nuclear dis- 
order strength is currently under evaluation as a potential early dia- 
gnostic marker of metastasis'^ and thus may be an indicator of early 
events in carcinogenesis. Measured values for have been found to 
be elevated with higher neoplastic behavior in genetically altered but 
microscopically indistinguishable human colon cancer cell lines 
One benefit of this kind of trans-network study is the ability to look 
at relationships between observations. For example, the AFM data 
presented suggests that metastatic cells are softer. Interestingly, the 
TFM data suggests that metastatic cells are able to exert increased 
force. This compound observation is both phenotypically intuitive 
and surprising. 

One of the unanticipated findings was that motility of metastatic 
MDA-MB-231 cells is profoundly different than that of MCF-lOA 
cells. MDA-MB-231 cells, in spite of exhibiting a slower overall 
speed, traveled farther and more linearly in the more pathologically 
relevant 3D motility assay than MCF-lOA cells, which move in a 
circular manner. Previously, motility measurements for these cells 
were performed in a 2D environment; in this milieu, MCF-lOA cells 
traveled farther than MDA-MB-231 cells. Motility was further 
explored through measurements of adhesion and rolling on matrices 



with different surface ligands. The MDA-MB-231 cell surface is very 
different from that of MCF-lOA cells, and these differences impact 
interactions with the microenvironment. CD44 distribution on the 
two cell lines differs; MCF-lOA cells exhibit abundant CD44 at the 
cell-substrate interface, whereas the MDA-MB-231 cells had less 
detectable CD44 at the cell- substrate interface. In addition, lipid rafts 
were detectable on the metastatic cells but not on the MCF-lOA cells. 
These results suggest that although the general topography of the two 
cell types is similar, the surface presentation of the cell adhesion 
protein CD44 and formation of lipid raft domains containing addi- 
tional cell surface receptors differ. These results are consistent with 
previous findings of elevated levels of cholesterol-rich lipid rafts in 
breast and prostate cancer cells compared with their normal counter- 

The metastatic cells exhibited less adhesion and rolling on E-selec- 
tin and HA surfaces, consistent with the lower expression of CD44 
detected at the substrate/membrane contact areas in TIRF images of 
the cells. The labeling experiments support the idea that the general 
topographies of the two cell lines are similar. The exclusion of CD44 
from the near-surface contact area in MDA-MB-231 cells or lack of 
proper post-translational CD44 modification may explain their lack 
of binding to the HA micro-patterned surface. These results differed 
substantially from those of Afify et aV^ who found that CD44 was 
required for breast cancer cell adhesion since anti-CD44 antibodies 
disrupted cell adhesion to HA-coated wells. In their study, Afify et al. 
only investigated the adhesion of two breast cancer cells lines: an 
MDA-MB-468 cell line that expresses high levels of CD44, and 
a T47D cell line with low to negligible expression of CD44; 



Table 2 Comparative experimental physical parameters for non-tumorigenic and metastatic cells 


Cell Physical Parameter 


Non-tumorigenic MCF-lOA 


Metastatic MDA-MB-231 


Morphology 


Cell shape and volume 

[2D°, 3D, and nuclear shape)^ 
Nuclear architecture 

Effect of matrix stiffness on cell proliferation 

and morphology 
Cell surface (CD44 expression patterns and 

lipid raft distribution) 


• Smaller cell volume, spherical cell shape 

• Heterogeneous, non-spherical nuclear shape 
Baseline 

Sensitive (increased proliferation with 
increased stiffness) 

• Abundant CD44 at cell surface-substrate 

interface 

• No measurable lipid rafts 


• Larger cell width/length ratio and volume 

• Kidney bean, non-spherical nuclear shape 
80% increased nuclear disorder 
Insensitive 

• Reduced CD44 at cell surface-substrate interface 

• Detectable lipid rafts 


Motility 


1 D motility in fibronectin channel 
2D motility on collagen 
3D motility in collagen 
Traction force on laminin 
Rolling on E-selectin 
Adhesion on hyaluronic acid 


Less displacement 

More displacement 

Less displacement (circular motion) 

Less traction force 

Yes 

Yes 


More displacement 

Less displacement 

More displacement (linear motion) 

More traction force 

No 

No 



Mechanics 



Deformability, elasticity (AFM Less deformable, less elastic More deformable, more elastic 

nanoindentation) 

Internal fluidity (microparticle displacement) Less fluidity More fluidity 

Stress response and survival 



Hypoxia 



Low pH environment 
Paclitaxel 



• Low cell viability, slow population recovery 

• No significant change in O2 consumption 

• Low nominal CEA expression v^ith small 

increase 
Greater toxicity 
Greater sensitivity 



• Low cell viability, subpopulation survival with fast 

recovery and greater proliferation 

• Fourfold reduction in O2 consumption 

• Low nominal CEA expression with moderate 
increase 

Sensitive 
Sensitive 



"Abbreviations used: 1 D: one-dimensional; 2D: 2-dimensionional; 3D: 3-dinnensionional; CEA: carcinoembryonic antigen. 
'^See text for additional explanation of experimental results. 
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MDA-MB-23 1 cells were shown to express CD44, but at a lower level, 
and their adhesion was not examined. In the Afify et at. analysis, HA 
was presented singularly, whereas in our studies HA was presented 
adjacent to PEG-ylated regions. Some of the present authors have 
previously shown that adhesion preference to HA-presenting region 
depends on the presentation of adjacent molecules^°. 

When subjected to stressors such as low O2 and acidosis, both 
MDA-MB-23 1 and MCF-lOA cells exhibited reduced viability as 
compared with culture under ambient O2 and physiological pH; 
MDA-MB-23 1 cell numbers were maintained to a greater degree 
than MCF-lOA, illustrative of the metabolic flexibility exhibited by 
metastatic cells compared to the non-tumorigenic cell line. MCF- 
lOA cells tolerated low O2 better in 2D as compared with 3D culture, 
with total cell numbers dropping by more than 50% after six days 
of culture in 3D. This finding was unexpected, and may point 
to different mechanisms of stress response mediated by culture 
dimensionality. In hypoxic conditions, the four-fold reduction in 
O2 consumption by the 3D-cultured MDA-MB-23 1 cells, in contrast 
to the unchanged MCF-lOA consumption, may be due either to 
increased phenotypic plasticity or to the preexistence of an emergent 
subpopulation with reduced O2 metabolism that survive under hyp- 
oxic stress. These results demonstrate that it is important to monitor 
not only the whole cell population but also single cells. The single cell 
studies of oxygen consumption support the idea that a subset of 
MDA-MB-23 1 cells restored the population after these microenvir- 
onmental stress-caused evolutionary bottlenecksa hypothesis that is 
also consistent with the results of the viability studies. Taken 
together, these results confirm and extend prior observations of 
the inherent capacity of metastatic cells to respond to microenviron- 
ment stressors of low 02 partial pressure and acidic pH. 

Proteomic network analysis of MCF-lOA and MDA-MB-23 1 cells 
confirmed changes in cell morphology and motility and identified 
genes that are involved in multiple metastatic hallmark features. For 
example, LOX, which is regulated in the morphology and survival 
networks by FBNl and HIFloc, respectively, is involved in extracellular 
matrix remodeling and is a hypoxia-responsive gene that represses E- 
cadherin, leading to cellular transformation and invasion. Through a 
combination of regulatory network inferences and proteomic studies, 
we interrogated the molecular origins of several classic phenotypic 
"hallmark" properties of tumor progression and metastasis. 

Clearly, an important next step is to bring the biophysical assays 
presented in this paper to the in vivo realm. Some of these assays 
could be modified and readily extended to the in vivo case, such as 
quantitative single-cell migration, traction force microscopy, and 
particle tracking microrheology. Other measurements cannot be 
readily conducted in vivo, such as 3D cytometry. 

In sum, the results of experiments comparing parameters of cell 
lines chosen to model the non-tumorigenic and metastatic states 
generated internally consistent results. The metastatic cells were 
more physically flexible and viable under a wider range of stresses 
than the non-tumorigenic cells. In addition, cell surface protein 
expression was related to the motility parameters measured. This 
pilot study validates the laboratory network approach and the use 
of physical sciences techniques to investigate cancer cell biology. 

Methods 

Cell lines. Cells were grown with standardized media and culture conditions 
described by Guise et al.^^ for MCF-lOA cells and Debnath et al.^ for MDA-MB-23 1 
cells. Additional specific cell growth conditions for specific assays are provided in the 
supplementary information (SI Methods). 

Cell morphology assays. Specific assay conditions are provided in SI Methods. 

Cell motility and mechanics assays. Specific assay conditions are provided in SI 
Methods. 

Cell stress response and survival assays. Specific assay conditions are provided in SI 
Methods. 



Gene and protein expression and analysis. Specific assay conditions are provided in 
SI Methods. 
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